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ABSTRACT 

Infrared ^^5-35 /im spectra for 40 solar-mass T Tauri stars and 7 intermediate-mass Herbig Ae stars with 
circumstellar disks were obtained using the Spitzer Space Telescope as part of the c2d IRS survey. This 
work complements prior spectroscopic studies of silicate infrared emission from disks, which were focused 
on intermediate-mass stars, with observations of solar-mass stars limited primarily to the 10 /im region. The 
observed 10 and 20 //m silicate feature strengths/shapes are consistent with source-to-source variations in grain 
size. A large fraction of the features are weak and flat, consistent with /im-sized grains indicating fast grain 
growth (from 0.1-1.0 /im in radius). In addition, approximately half of the T Tauri star spectra show crystalline 
silicate features near 28 and 33 /im indicating significant processing when compared to interstellar grains. A 
few sources show large lO-to-20 /im ratios and require even larger grains emitting at 20 /im than at 10 /im. 
This size difference may arise from the difference in the depth into the disk probed by the two silicate emission 
bands in disks where dust settling has occurred. The 10 /im feature strength vs. shape trend is not correlated 
with age or Ha equivalent width, suggesting that some amount of turbulent mixing and regeneration of small 
grains is occurring. The strength vs. shape trend is related to spectral type, however, with M stars showing 
significantly flatter 10 /im features (larger grain sizes) than A/B stars. The connection between spectral type 
and grain size is interpreted in terms of the variation in the silicate emission radius as a function of stellar 
luminosity, but could also be indicative of other spectral-type dependent factors (e.g. X-rays, UV radiation, 
stellar/disk winds, etc.). 

Subject headings: circumstellar matter — stars:pre-main sequence-infrared: ISM: lines and bands — stars: 
formation — solar system: formation 



1. EVITRODUCTION 

Dust in disks can be quite different from dust in the in- 
terstellar medium (ISM). Observations and chemical model- 
ing (Grossman 1972, Gail 1998) suggest that the dust in the 
early stages of star formation is primarily composed of small 
(<1 /im) amorphous silicates with strong features at approx- 
imately 9.7 and 18.5 /im (for a summary, see Pollack 1984). 
Large modifications of the dust occur in the envelopes and 
disks around young stars, as the initially small grains are pro- 
cessed via collisions and coagulation. Spectral energy dis- 
tributions (SED) indicate that grain growth, and the corre- 
sponding settling of large grains to the disk midplane, is oc- 
curring in some disks. Dust settling and growth affect disk 
temperatures and vertical structures, resulting in dust photo- 
spheres that are fl atter rat her than flared (e.g., D'Alessi o et alJ 
ll999HChiang et aL..200klDullemond & Dominik.2004 . Ad- 
ditionally, some main sequence stars show evidence of a "sec- 
ond generation" of small grains in a debris disk, produced by 
the collision and fragmentation of planetesimals, that may be 
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quite different from the "primordial" dust. 

The spectroscopic study of silicate emission has proven 
a valuable tracer of grain processing within young circum- 
stellar disks. This method probes small grains via optically 
thin emission from the surface layer of generally optically 
thick disks. Studies of ^10 /im silicate emission from Her- 
big Ae/Be stars (hereafter HAEBE) and T Tauri stars (here- 
after TTs) with disks show early evidence of the growth of 
irains from ~0.1 to 2.0 /im (see e.g.. 



these surface-l ayer g i 
'Bouwman et al.!'200r; 



'van Boekel et al.'*2003'). Detections of 
additional spectral features arising from crystalline silicate 
emissi on in som e HAEBE and TTs disks ( Waelkens et aD 
J996; Sitko et al. 1999; Meeus et alJ200UlHonda et aL2003F 
Ucke & van den Ancker 2004), the debris disk Pictoris 
(e.g., Knacke et al. 1993), and comets (see review by 'Wooden] 
1^02) provide evidence for silicate processing during the disk 
phase. Although changes in grain size and composition are 
closely linked to disk properties and planet formation, the rate 
and mechanism of grain growth and processing in disks are 
still not well understood. 

Previous studies of silicate emission from TTs disks were 
primarily focused on the Si-0 stretching mode feature near 
10 /im, which can be observed from the ground. Grain growth 
and crystallization, however, have similar effects on the shape 
of the 10 /im feature. In these previous studies, the pres- 
ence of crystalline silicates was often established through 
the presence of flux at 11.3 /im, corresponding to an emis- 
sion feature of the crystalline Mg-rich silicate forsterite, in 
addition to the presence of an amorphous olivine feature at 
9.8 /im. The combined effect is a broad, trapezoidal sili- 
cate emission feature with a peak near 9.7 /im and a sec- 
ondary peak near 11.3 /im. If crystalline silicates are abun- 
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dant, the feature can appear broader and flatter. In con- 
trast, growth of amorphous oHvines can mimic the effects of 
crystafline-amorphous silicate mixtures, resulting in weaker, 
"flat-topped" 10 nm silicate features, with s imilar flux at 9.8 
and 11.3 (e.g., jv an Boekel et al. 2003; Przvgodda et ajj 
l200 3t iKessler-Silacci et al. 2005). Indeed, the above studies 
showed that the flattened silicate features were well matched 
by models of purely amorphous olivines with grain sizes 
(radii) of 2 /im. Therefore, evaluation of the grain size and/or 
crystallinity requires observations of a larger spectral region, 
in which the presence of distinct forsterite or enstatite features 
(and PAHs, which emit near 7.7, 8.5, and 1 1.2 /im) can be as- 
sessed. Complementary studies of the O-Si-O bending mode 
near 20 /im probe slightly cooler dust and combined give a 
better sense of the grain size and crystalline fraction in the 
disk. 

Spectra covering a large spectral region, including sev- 
eral isolated crystalline silicate features (27.5, 33.5, 35.8 
and 70 /im), were obtained for disks around HAEBE 
stars (and some TTs) with the Infrared Space Obser- 
vatory (ISO). ISO observations of several disks around 
intermediate-mass HAEBE stars clearly esta bhshed the pres- 
ence/absence of crystalline silicates (see iMolster & WatersI 
[2003; vanDishoeck 2004, for reviews). Additionally, ISO 
observations probed the ~10 and ~20 emission features 
from amorphous silicates in HAEBE disks, allowing grain 
growth to be studied (see Acke & van den Ancker 2004, for 
a review). Due to sensitivity limitations, ISO studies of sil- 
icates focused primarily on intermediate or high-mass stars. 
Studies of TTs have therefore been limited, for the most part, 
to the 10 /im region. Although changes in the shape of the 
10 /im em ission features are observed in small sampl es of 
TTs ( Przvgodda et aljl2003t IK^^-Silacci et al.ll200l. the 
cause of these changes is unclear and could be a combination 
of silicate grain growth and crystallization. 

In this study, the improved sensitivity of the Infrared Spec- 
trograph (IRS) aboard the Spitzer Space Telescope was used 
to expand such spectroscopic studies to include a large sam- 
ple of disks around low-mass, sun-like stars, creating a 
database analogous to ISO studies of '^high/intermediate- 
mass stars. A few HAEBE stars were also observed for com- 
parison. The data presented here are part of the c2d Spitzer 
legacy program designed to study the evolution of circumstel- 
lar matter 'From Molecular Cores to Planet-Forming Disks' 
llEvansetalJl200 3). This paper contains preliminary results 
from the program and will focus on statistical analysis of the 
strongest silicate features. It will be followed by a more de- 
tailed study of weaker features, including emission from crys- 
talline silicates. The Spitzer IRS observations are described in 
§2. In §3, an inventory of the most prominent silicate emis- 
sion features will be presented. The observed spectra will be 
compared with models of amorphous silicates in §4. We per- 
form a statistical analysis of the 10 and 20 /im silicate fea- 
tures and interpret these trends in terms of grain growth in §5. 
Then in §6, we examine the relationship between the strength- 
shape trends found in §5 and spectral type, stellar age, and Ha 
equivalent width. Finally, in §7, we discuss the crystalline sil- 
icate emission features in the observed spectra. 

2. OBSERVATIONS 

Spectra for 40 solar-mass T Tauri stars and 7 intermediate- 
mass HAEBE stars were obtained using the Infrared Spec- 
trograph (IRS) aboard the Spitzer Space Telescope from De- 
cember 2003 through December 2004. The observations and 



source parameters are described in Table ^ for the T Tauri 
stars (top) and HAEBE stars (bottom) in our sample. The 
spectra were observed with combinations of the short-low 
(SL), short-high (SH), long-low (LL), and/or long-high (LH) 
modules (see Tablencolumn 4). SL (A — 5.3-14.5 /im) and 
LL (A — 14.2^0.0 /im) spectra have a resolving power of 
R = A/AA - 100 and the SH (A = 10.0-19.5 /xm) and 
LH (A ~ 19.3-37.0 /im) spectra have a resolving power of 
^600. For approximately half of the sources, SL spectra are 
part of the GTO programs and not yet available for analysis 
and, therefore, there is no spectral information short-ward of 
10.0 /im. Exposure times were chosen to achieve signal-to- 
noise ratios of 50 and 100 for sources brighter and fainter than 
500 mJy, respectively, with the weakest sources in our sam- 
ple having fluxes of 100-200 mJy at 15 /im. All sources in 
the c2d program (Evans et al. 2003) that were observed prior 
to December 2004 and show evidence of silicate emission are 
included (47 sources; see TableQ. 

Data reduction was done via the c2d Interactive Analy- 
sis (c2dia) reduction environment. C2dia contains optimized 
extraction algorithms developed by the c2d legacy team^. 

For wavelength calibration and IRS apertu re definition, 
tools from the SMART'" software package (iHigdon et alJ 
120041) were used. The extracted spectra were defringed 
using the IRSFRINGE" package developed by the c2d 
team (Lahuis & Boogert 2003). The spectra were extracted 
from the Spitzer Science Center (SSC) Basic Calibrated Data 
(BCD), pipeline version SI 1.0.2'^. Two different extraction 
methods were used and compared to reduce spectral artifacts 
resulting from the extraction process. 

The first method uses a full aperture extraction for SH and 
LH and fixed width aperture extraction for SL and LL. The 
extraction aperture was chosen to be large enough to enclose 
the complete source. (At the short wavelength end, the ex- 
traction aperture is wider than that used in the SSC pipeline. 
However, since the SL and LL spectra are corrected for the 
background, as described below, this has a negligable effect 
on the extracted spectrum.) Bad pixels are corrected by in- 
terpolating in the cross-dispersion direction using a fit to the 
order-averaged source profile. For SL and LL the the large 
apertures allow off-source spectra to be extracted and these 
spectra were used to make wavelength dependent background 
corrections. No background data were available for SH and 
LH spectra. A spectrum was extracted for each position, com- 
pared to check for artifacts, and finally averaged to produce 
the final spectrum. 

In the second method, all BCDs from both dither positions 
are combined within the extraction algorithm. The extrac- 
tion is performed by integrating over a source profile fit in 
the cross-dispersion direction. The source profile is a tem- 
plate created from standard star observations (including sky 
measurements for SH and LH). The width and center of this 
template is adjusted for each observed source in our sample to 
encompass 95% of the observed flux. This fitting process is 

' The c2d extraction algorithms will become publicly available through the 
Spitzer Space Science Center as part of the c2d legacy program. 

The SMART software package is publicly available on the website 
\protect ittp://ssc. spitzer.caltech.edu/archanaly/contributed/smart 

" IRSFRINGE is included in SMART 

and available as a stand-alone package from 
IXprote ct ittp://ssc. spitzer.caltech.edu/archanaly/contributed/irsfringe/ 

'■^ Recently data from the SSC pipeline version S12 as well as more ad- 
vanced versions of the extraction algorithms have become available. We have 
verified that these do not have a significant impact on the reduced spectra and 
do not affect the results presented in this paper 
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TABLE 1 
Source list 





RA 


DEC 


wDserved 


AOR 


UDserv3.tion 


Age 


nuLCj VV 


opectrsl 




Source 


(JZUUU) 


(J zUUU J 


Modules 


J\.ey 


Date 


(Myr) 


(A) 


Type 


alt 

Kei. 


RNO 15 


03 27 47.68 


-f 30 12 04.3 


SL SH LH 


0005633280 


2004-08-30 




116 




1 


LkHa 327 


03 33 30.41 


+31 10 50.4 


SL SH LH 


0005634560 


2004-02-03 


1.4 


51-65 


K2 


1,2,3 


LkHa 330 


03 45 48.29 


+32 24 11.8 


SL SH LH LLl 


0005634816 


2004-10-25 




11-20 


G3 


2,3 


IRAS 03446+3254 S 


03 47 47.12 


+33 04 03.4 


SL SH LH LLl 


0005635072 


2004-09-29 










VVlOTau'' 


04 31 57.79 


+ 18 21 36.3 


SH LH 


0005636608 


2004-09-29 


0.7/0.4 


11/89 


M0.5/M3 


4 


CoKu Tau 4 


04 41 16.79 


+28 40 00.5 


SH LH 


0005637888 


2004-09-02 


1.2-1.6 


1.8-2.8 


M1.5 


2,5,6 


IRAS 08267-3336 


08 28 40.70 


-33 46 22.3 


SL SH LH LLl 


0005639168 


2004-11-11 




25-35 


K2-K3 


7,8 


SX Cha 


10 55 59.74 


-77 24 39.9 


SH LH 


0005639424 


2004-08-31 


1^ 


26.7 


M0.5 


9-11 


SY Cha 


10 56 30.47 


-77 1 1 39.4 


SH LH 


0005639424 


2004-08-31 


2-5 


24-64 


MO 


9-12 


TWCha 


10 5901.11 


-77 22 40.8 


SH LH 


0005639680 


2004-09-01 


20 


26.1 


MO 


9,11 


VW Cha'' 


11 08 00.53 


-77 42 28.7 


SH LH 


0005639680 


2004-09-01 


0.4-0.9 


72-147 


K2 


9-12 


VZCha 


11 09 23.80 


-76 23 20.7 


SH LH 


0005640448 


2004-09-02 


4-10 


58-71 


K6 


9-12 


WX Cha'' 


11 09 58.75 


-77 37 08.9 


SH LH 


0005640192 


2004-09-01 


1-5 


65.5 


K7-M0 


9-11 


ISO Cha237 


11 10 11.44 


-76 35 29.2 


SH LH 


0005640448 


2004-09-02 




<3 


MO 


13 


C7-11'' 


11 10 38.01 


—77 32 39.9 


SH LH 


0005640192 


2004-09-01 


0.2-1 


4.0 


K3 


9,14 


HM 27 


11 10 49.62 


—77 17 51.7 


SH LH 


0005640192 


2004-09-01 


30-40 


200.0 


K7 


9,11 


XX Cha 


1111 39.67 


—76 20 15.1 


SH LH 


0005640448 


2004-09-02 


2^0 


133.5 


M2 


9-^11 


T Cha 


11 57 13.53 


-79 21 31.5 


SH LH 


0005641216 


2004-07-18 


>12.5 


2-10.0 


G2-G8 


15-17 


IRAS 12535-7623 


12 57 11.78 


-76 40 11.5 


SH LH 


0011827456 


2004-08-31 




3.0 


MO 


17,18 


Sz50 


13 00 55.37 


-77 10 22.2 


SH LH 


0011827456 


2004-03-25 




46 


M3 


19 


HT Lup'' 


15 45 12.87 


-34 17 30.6 


SL SH LH 


0009829120 


2004-08-28 


0.4-0.8 


3-7 


K2 


11,12,20 


GWLup 


15 46 44.68 


-34 30 35.4 


SL SH LH LLl 


0005643520 


2004-08-30 


1.3-3.2 


90-98 


M2-M4 


11,20,21 


HMLup 


15 47 50.63 


-35 28 35.4 


SLLLl LL2 


0005643776 


2004-08-30 


1.2-2.7 


115-155 


M3 


11,20 


Sz 73 


15 47 56.98 


-35 14 35.1 


SL SH LH LLl 


0005644032 


2004-08-30 


2.6-5.4 


97-150 


K2-M0 


11,20,21 


GQLup 


15 49 12.10 


-35 39 05.0 


SL SH LH LLl 


0005644032 


2004-08-30 


0.1-0.6 


31-39 


K7-M0 


11,12,20 


IMLup 


15 56 09.17 


-37 56 06.4 


SL SH LH LLl 


0005644800 


2004-08-30 


0.1-0.6 


4.7-8.1 


MO 


11,12,20 


RULup 


15 56 42.31 


-37 49 15.5 


SL SH LH LLl 


0005644800 


2004-08-30 


0.1-0.5 


136-216 


K7-M0 


11,12,20 


RYLup 


15 59 28.39 


-40 21 51.2 


SL SH LH LLl 


0005644544 


2004-08-30 


1.6-3.2 


7.3 


K0-K4 


12,20 


EX Lup 


16 03 05.52 


-40 18 24.9 


SL SH LH LLl 


0005645056 


2004-08-30 


1.4-3.0 


31^3 


MO 


11,12,20 


Sz 102 


16 08 29.70 


-39 03 11.3 


SL SH LH 


0009407488 


2004-03-25 




377 


K0-M4 


20 


AS 205'' 


16 11 31.35 


-18 38 26.1 


SL SH LH 


0005646080 


2004-08-28 


0.1 


55-155 


K5/M3 


2,3,12,22 


VSSGl 


16 26 18.86 


-24 28 19.7 


SH LH 


0005647616 


2004-08-28 










DoAr 24E'' 


16 26 23.38 


-24 21 00.1 


SH LH 


0005647616 


2003-12-15 


1.5 


5 


KO 


22,23 


GY23 


16 26 24.06 


-24 24 48.1 


SH LH 


0005647616 


2004-08-28 






K5-M2 


24 


SR 21 N 


16 27 10.28 


-24 19 12.5 


SH LH 


0005647616 


2003-12-15 


1 


0.54 


G1-G2.5 


22,25 


SR 9^ 


1 6 27 40 27 


—24 22 04.0 


SH LH 


001 2027392 


2004-09-02 




6-14 


K5-M2 


1 23 24 26 


Haro 1-17 


16 32 21.94 


—24 42 14.7 


SL SH LH LLl 


001 1827712 


2004-08-29 




15 


M2.5 


26 


RNO 90 


16 34 09.18 


-15 48 16.8 


SL SH LH LLl 


0005650432 


2004-08-28 


6 


76 


G5 


1 


EC 82 


18 29 56.89 


+01 14 46.5 


SL SH LH 


0009407232 


2004-03-27 




5-11 


MO 


2,27 


CK4 


18 29 58.21 


+01 15 21.7 


SL SH LH 


0009407232 


2004-03-27 










BFOri 


05 37 13.26 


-06 35 00.6 


SL SH LHLLl 


0005638144 


2004-10-03 


2.0 


6-11 


A5-F6 


2,12,28-30 


RRTau 


05 39 30.52 


+26 22 27.0 


SL SH LHLLl 


0005638400 


2004-09-28 


0.6 


21.2-50 


B8-A5 


2,28,29,31 


HD 98922 


11 22 31.67 


-53 22 11.4 


SHLH 


0005640704 


2004-01-04 




27.9 


B9 


30,31 


DLCha 


13 06 08.36 


-77 06 27.3 


SHLH 


0005642240 


2004-07-14 






M6 


22 


HD 135344 


15 15 48.44 


-37 09 16.0 


SHLH 


0005657088 


2004-08-08 


17 


17.4 


A0-F4 


30,32-34 


HD 163296 


17 56 21.29 


-21 57 21.9 


SHLH 


0005650944 


2004-08-28 


4-6 


12-17 


A0-A2 


15,29-32 


VV Ser 


18 28 47.86 


+00 08 39.8 


SL SH LH 


0005651200 


2004-09-01 


0.6 


22-90 


B1-A3 


3,28-31 



References. — (I) Levreault 1988; (2) Cohen & Kuhi 1979; (3) Fernandez et al. 1995; (4) Hartiean et al. 1994; (5) D'Alessio et al. 2005; (6) Kenvon et al. 1998; (7) 
iReiminh & Pettersson 1993; (8) Sahu & Sahu 1992; (9) Lawson et al. 199fi; (10) Hartmann et al. 1998; (11) Ap nenzelle r et al. 1983; (12) Reiouith et al. 1996; (13) Saffe etall 
ISXH: (14) Hartiean 1993; (15) van den Ancker et al. 1998; (16) Alcala et al. 1995; (17) Gregorio-Hetem e t alil992t7 l8) Saitori et al. 2003; (19) Hughes & Hartigan 1223 
(20)'HuKhes et al. 1994; (21) Heyer & GrahaiTi 1 989; (22) P ra j o e t ah 2003 ; (23) Bouvier & ADpenzelle iil99a (24)jLuh man & Rieke 1999; (25) Marti n et al. 1998 - (26) 
iRvdsren 1980; (27) Gomez de Castrcll99^ (28) lNatta et alll99Tt (29) lFi'n'kenzeller & Mundtl984t (30) IAcke et all2003t (31) lThe et aHl994 (32) lThi et ail200lt (33) 
[Houk 1982; (34) Dunkin et al. 19971 

Note. — Sources in the top portion of the table are T Tauii stars; sources in the bottom portion are HAEBE stars. 

^AU Ha lines are in emission, except that of SR 21, which is in absorption 

''This source is a binary in 2MASS K-band images, with a separation < 3", and is unresolved with Spitzer IRS. Spitzer IRS observations are pointed at the center of the 
binary system. Stellar/disk parameters quoted from the literature include both sources 



performed for each source by using the highest quahty data. 
Once the source profile is fit, a (uniform) local sky back- 
ground level, which is wavelength dependent, is estimated. 
This method also reduces the effects of unidentified bad pix- 
els. In particular, in cases where the BCD images are largely 
affected by bad pixels (e.g., for LH), this method quite often 
gave significant improvements over the full aperture extrac- 



tion. In the version of c2dia used here, the wavelength and 
flux calibration files produced by SSC pipeline were used. For 
the S12 version of c2dia, the extraction algorithm is further 
developed and independently calibrated using a suite of stan- 
dard star calibrators. The extraction algorithms and the cali- 
bration involved will be described in more detail in Lahuis (in 
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prep.). 

To correct for the possibly significant sky contribution in 
the SH and LH spectra, spectra obtained using full aperture 
extraction were corrected using the background estimate from 
the source profile fitting extraction method. Both extractions 
were also compared against the spectra extracted from the 
SSC pipeline. In a few cases, particular modules suffered 
from artifacts due to data selection and/or bad-pixel correc- 
tion and the corresponding modules of the SSC spectrum were 
used. In cases where the SSC spectra were used, a correction 
to the zero level was applied to correct for background emis- 
sion. Good sections of the background corrected spectra were 
used to estimate the zero level in the other modules. When 
necessary, individual orders and/or modules were multiplied 
by small (< 5% and < 15%, respectively) constant factors to 
correct flux offsets between orders/modules. These flux off- 
sets are likely related to pointing errors resulting in the source 
not being centered in the aperture and seem to be corrected in 
the SSC pipeline version S12. 

3. INVENTORY OF SILICATE EMISSION FEATURES 

The extracted spectra of TTs and HAEBE disks are shown 
in Figure [0 and Figure |2l respectively. Two sources, HT Lup 
and HD 163296, were observed twice, once in IRS-staring 
mode and once in spectral mapping mode. The two observa- 
tions for each source do not differ significantly; integrated in- 
tensities, peak-to-continuum fluxes, and feature widths agree 
to within ^1% for the 10 and 20 /im features. Therefore, only 
one spectrum for each source is included in Figures [2 and |2] 
and Table 121 

An inventory of the silicate emission features is presented in 
Table|2] The features listed are relatively isolated and can be 
clearly identified in the spectra. The peak wavelengths of the 
silicate features depend on the exact composition and grain 
size of the dust. Therefore, the broad amorphous silicate fea- 
tures, corresponding to the Si-O stretching and O-Si-O bend- 
ing modes, are labeled as '^10 and ~20 /im. The peaks of 
the identified crystalline silicate features may vary within the 
wavelength ranges given in Table |2] which correspond to the 
ranges in peak wavelengths for features of crystalline enstatite 
and forsterite of sizes between 1 and 10 /im. Features are 
marked with "Y", "T", or "N" if the feature is detected, tenta- 
tively detected, or absent. 

Emission at 10 and 20 /im from amorphous silicates dom- 
inates the spectra of most of the disks observed. Broad 
20 yum amorphous silicate features can be seen in almost 
all of the spectra, while the 10 /im amorphous silicate fea- 
tures are slightly less prominent (seen in 42 of 47 spectra). 
Emission near 11.3 /im, due to PAH or crystalline forsterite, 
is clearly identified in ^1/3 of the spectra. This is a lower 
limit, however, as tentative emission is seen in several addi- 
tional spectra in which the presence of artifacts prevent fur- 
ther analysis. IS O observa tions of disks around HAEBE stars 
(^cke & van de n Anckeil2 004) report 11.3 /im emission from 
~l/4 of the observed sources, but the same studies find that 
the fraction showing PAH emission at 7.7 /im is much higher 
(57%). 

The lattice modes of crystalline silicates at A > 25 /im ap- 
pear to be more prominent than the 11.3 /im feature in our 
data and do not suffer from confusion with PAH emission. 
Approximately half of the disks in our sample show emission 
features near 28 /im and in the 33-35 /im region, characteristic 
of crystalline enstatite and forsterite. The 24 to 36 /tm region 
is shown for 6 of these sources in Figure |3j with shaded re- 



gions depicting the typical positions of enstatite and forsterite 
features. These data indicate that crystalline silicates can be 
quite prominent in TTs disks. This is in surprising contrast to 
ground-based studies of silicates in the 10 /im region, which 
found crystalline silicate emission from very few TTs (see 
.Kessler-Silacci et al. 2005 , for summary). 

As the 10 and 20 /tm emission features from amorphous 
silicates are the most easily identified and isolated, the rest of 
this paper will focus on the analysis of these features. PAHs 
will be discussed in detail in Geers et al. (in prep.) and SEDs 
will be modeled in Augereau et al. (in prep.). Silicate com- 
position will be modeled in more detail by Kessler-Silacci et 
al. (in prep.), including analysis of all crystalline silicate fea- 
tures. 

4. COMPARISONS WITH MODELS OF AMORPHOUS SILICATES 

As mentioned above, the peak wavelengths, strengths and 
shapes of the silicate features depend on the exact composi- 
tion, sizes, and shapes of the silicates. Therefore, in order to 
interpret the observed silicate emission, we must compare our 
data with models of the emission from a variety of silicate 
grains. Modeled grain opacities can be compared to the ob- 
served emission once the continuum is removed. In §4.1, we 
describe the method of continuum normalization used here. 
Our procedure for modeling of grain opacities from labora- 
tory optical constants is described in §4.2. Finally, in §4.3, we 
will quantitatively compare modeled and observed opacities 
in the 10 and 20 /tm regions. 

4.1. Continuum normalization 

Continuum normalization was performed using one of three 
methods, as noted in Table|2] and described as: 

1. For spectra that include SL (e.g., RNO 15), a 2nd order 
polynomial was fit to the continuum in the following 
regions: 6.8-7.5 /im, 12.5-13.5 /im, and 30-36 /tm. 

2. For spectra that do not include SL (e.g., V710 Tau), the 
continuum was obtained with a linear fit to the 12.5- 
13.5 /im and 30-36 /im regions. 

3. In cases where the SED appears to fall sharply with in- 
creasing wavelength short-ward of 12-15 /tm and then 
rises sharply beyond '^15 /tm (e.g., LkHa 330), we fit 
the 10 and 20 /im features separately, using a 2nd order 
polynomial fit for wavelengths short- ward of ~15 /im 
and a linear fit long-ward of ^15 /im. These sources 
will be discussed in more detail in Brown et al. (in 
prep.). 

The continuum fits are overlaid on the observed spectra in 
Figures ^ and |2] as dash-dotted lines. The continuum fit- 
ting regions were chosen to represent the areas least affected 
by features in the spectra and, as shown, vary slightly from 
source to source. Opacities obtained from laboratory spec- 
tra of amorphous olivines (see §4) indicate that the regions 
chosen above should be clear of amorphous silicate emission 
features, with the exception of the region from ~13-15 /im, 
which likely contains contributions from the wings of both 
the 10 and 20 /im features. Thus, the feature strengths ob- 
tained from the continuum normalized spectra will be under- 
estimated. 

After the continuum fitting was performed, the normalized 
spectra 5^ were obtained via application of the formula. 
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Fig. 1 . — Spitzer mid- to far-IR spectra for T-Tauri disks. Tlie spectra for eacli source are shown in units of Log XFx. The SH and LH spectra have been 
median smoothed over 3 and 5 channels, respectively. No smoothing has been applied to the other modules. The dash-dot line denotes the fit to the continuum 
described in the text and Tablel2l 
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Fig. 2. — Spitzer mid to far-IR spectra for HAEBE disks. The spectra for each source are shown in units of Log XFx. The SH and LH spectra have been 
median smoothed over 3 and 5 channels, respectively. No smoothing has been applied to the other modules. The dash-dot line denotes the fit to the continuum 
described in the text and Tablel2l 



where Fi^ is the observed spectrum and i^^ ^ is the fitted con- 
tinuum, both in units of Jy, and < F^^c > is the frequency 
averaged continuum flux. We divide by a frequency aver- 
aged value of the continuum flux in order to remove any de- 
pendence on the slope of the continuum. This normaliza- 
tion method, in cluding the offset of 1 .0, is c onsistent with the 
method used bv ivan Boekel et all J2003I) and lPrzv^odda et aO 



("2003^, to which we compare in §5. As the continuum may 
change dramatically between 10 and 20 /im for some sources 
(e.g., CoKu Tau 4), < i^^.c > is calculated separately for re- 
gions of A = 5-13 /im and A = 13-37 /im. If c is used 
in place of < c >, then the normalization would reduce 
to F^/F,y^c, which is approximately equal to the optical depth 
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TABLE 2 
Silicate emission features 



Species 


Silicates" 


PAH/Cr.Sil. 


Silicates" 


Cr.Sil.'^ 


Cr.Sil.'' 


Norm.'= 


X(p.m) 


~10 


11.2-11.3 


~20 


28-29 


33-35 


Metliod 


RNO 15 


Y 


T 


T 


Y 


T 


\ 




Y 


Y 


Y 


Y 


Y 


J 




Y 


Y 


Y 


T 


isr 


J 




Y 


Y 


Y 


N 


T 


J 


V71 Tnii 


Y 


T 


Y 


T 


T 


2 


r'nPCii Tnii A 

^^UrvLI 1U.L1 '-T 


Y 


T 


Y 


N 


N 


2 




Y 


N 


Y 


N 


N 


J 




Y 


T 


Y 


T 


Y 


2 


SY Cha 


Y 


T 


Y 


T 


T 


2 


TW Cha 


Y 


T 


Y 


T 


T 


2 


VW Cha 


Y 


T 


Y 


Y 


Y 


2 


VZ Cha 


Y 


T 


T 


Y 


Y 


2 


VV A V^Ild 


Y 


T 


Y 


T 


Y 


2 


ISO Cha237 


Y 


T 


Y 


T 


T 


2 


C7-11 


Y 


T 


Y 


T 


T 


1 


KM 27 


Y 


N 


Y 


T 


T 


2 




Y 


N 


Y 


T 


T 


2 


T Cha 


N 


Y 


T 


N 


N 


3 




Y 


T 


Y 


Y 


N 


2 


Sz 50 


Y 


N 


Y 


N 


T 


2 


HT Lup 


Y 


T 


Y 


T 


Y 




nW T iin 


Y 


T 


Y 


T 


T 




nivi i^up 


Y 


T 


T 


T 


T 




Sz 73 


Y 


T 


Y 


T 


T 




no T iin 


Y 


T 


Y 


Y 


Y 




TM T iin 
liVl l_jUp 


Y 


T 


Y 


T 


T 




RTT T iin 


Y 


T 


Y 


Y 


Y 




RY Lup 


Y 


N 


Y 


N 


N 




FY T „Ti 


Y 


T 


Y 


T 


T 




Sz 102 


Y 


T 


Y 


N 


Y 




AS 205 


Y 


T 


Y 


T 


T 




VSSGl 


Y 


Y 


N 


Y 


T 


2 


DoAr 24E 


Y 


Y 


Y 


Y 


T 


2 


GY23 


N 


Y 


N 


Y 


T 


2 


SR21 


N 


Y 


Y 


T 


T 


3 


SR9 


Y 


T 


Y 


T 


N 


2 


Haro 1-17 


Y 


Y 


Y 


Y 


Y 


1 


RNO 90 


Y 


N 


Y 


T 


N 


1 


EC 82 


Y 


N 


Y 


N 


T 


1 


CK4 


Y 


T 


Y 


T 


Y 


1 



BFOri Y T Y T Y 1 

RRTau N Y Y Y Y 1 

HD 98922 T Y Y Y Y 2 

DLCha T Y Y Y N 2 

HD 135344 N Y T N N 3 

HD 163296 Y T Y T Y 2 

VVSer Y T Y T Y 1 



Note. — Features are marked witti a "Y" if detected, "N" if not detected, or "T" if the identification is 
tentative. 

"Broad, primarily amorphous olivine/pyroxene emission features. 

^Feature includes PAH at 11.2 /^m and crystalline forsterite at 11.3 ^m (Geers et al., in prep.). 
'^Feature includes crystalline forsterite near 27.9 /im and/or crystalline enstatite near 28.5 /jm 
^Feature includes crystalline forsterite near 33.6 /^m and/or crystalline enstatite near 34.5 ^m 
^Continuum normalization procedure as desciibed in the text 



for optically thin emission; 

Ij^^l-e-^'^ ^T, (r, «1), 

where is the optical depth as a function of frequency. Com- 
parison between the two continuum normalization methods 
for the spectra in the c2d sample shows that the 10 and 20 /im 
feature strengths and shapes do not differ significantly (~ 5% 
and ^ 10%, respectively). 



Although continuum normalization has been previ- 
ously used to interpret ground-based and I SO observa- 
tions of 10 /im silicate emission features ("Honda et alj 
20031 Ivan B oekel et al. 2003; Przvgodda et al. 200 J 
Acke & van den Anck er 2004), it was previously not applied 
to 20 /im siUcate features because the emission feature 
strength can be confused by a rising continuum. To address 
this issue, we plot the spectral index (a) of the derived 
continuum from 13 to 35 /im versus the strength of the 
continuum normalized 20 /im feature (5*^°^™) in Figure |3 
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Fig. 3. — Continuum normalized and smoothed LH spectra for 6 of the 
sources that show emission from crystalline silicates in the 26.5-30 and/or 
33-35 ^m regions (depicted by shading). Each pair of spectra contains one 
HAEBE star (bottom) and one T Tauri star (top). Features in the 26.5-30 fim 
region tend to be narrower for the T Tauri stars. 
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Fig. 4. — Evaluation of the continuum normalization method at 20 ^m. 
13-35 /im spectral indices are plotted versus the strengths of the continuum 
normalized 20 /im silicate emission features for all sources. TTs are denoted 
by filled squares and HAEBE are denoted by filled triangles. The mid-IR 
spectral indices and 20 fim feature strengths appear to be at most weakly 
inversely correlated (r = —0.2, 80% significance), with more steeply rising 
spectra corresponding to stronger 20 fim silicate emission features. 



The spectral indices are calculated as. 

Log (AfcFAj-Log {XaFxJ 



Ola-b 



Log (Afc) - Log (Aa) 



is the flux-density 
-1 The bulk of 



where Aa = 13 ^m, Af, — 35 /im, and Fx 
at wavelength A^; in units of erg cm^^ s^^ /im 
the spectra are relatively flat or falling in units of XF\, with 
a = 0-1 and 5^°;^" = 1.1-1.5 (see Figure gj. The five 
sources for which we use continuum fitting method 3 have 




c2d observations 
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Fig. 5. — Evidence of grain growth in the Si-O stretching and O-Si-O 
bending mode features. The top panels show the observed normalized spectra 
in the a) 10 fim and b) 20 fim regions for sub-samples of our sources. The 
bottom two panels show the normalized absorption efficiencies (Qabs) for 
models of spherical grains of amorphous olivines with various grain sizes 
calculated for the 10 fim and 20 ^m regions. Models of filled homogeneous 
spheres calculated using Mie theory are shown in c) and d) and models of 
hollow spheres calculated using DHS theory are shown in e) and /). Spectra 
in all panels have been artificially shifted along the y-axis by a constant value 
as a function of wavelength, such that the spectrum of each source could be 
seen more clearly. The minimum of each normaUzed spectra was 1 .0 prior to 
adding these constants. 



spectra that are steeply rising from 13 to 35 /im (a = —1 
to —2) and lie along the bottom of Figure |3 A few other 
sources with rising spectra (e.g., RY Lup, EC 82) can be seen 
in the right of Figure|3 We find that the spectral index and the 
derived strength of the 20 fim feature is at most weakly corre- 
lated (at the 80% significance level), with the more steeply 
rising spectra corresponding to the stronger 20 fim silicate 
features. There is thus a small chance that the shape/strength 
parameters derived for a particular source may be related to 
the continuum fit. As stated above, modeling of individual 
spectra without continuum normalization will be presented in 
Kessler-Silacci et al. (in prep.). 

A representative sample of continuum normalized spectra 
in the 10 and 20 fim regions are shown in Figure|5^,b. A large 
fraction of the observed spectra are similar to GQ Lup, pos- 
sessing weak, and sometimes flat-topped 10 fim silicate emis- 
sion features. A few others possess centrally peaked 10 fim 
features, with large peak-to-continuum flux, as is the case for 
EC 82. The rest of the sample falls between these two ex- 
tremes. Most spectra show 20 fim features with lower peak- 
to-continuum flux in the 20 fim region with respect to that at 
10 fim. 

4.2. Modeling dust optical properties 

To aid in the interpretation of the observed silicate emis- 
sion features, we model the opacities for a sample of grains 
of different shape, size and composition. Optical constants 
are obtained from laboratory data available from the Jena - 



8 



Kessler-Silacci et al. 



TABLE 3 
List of adopted grain material 



Dust type 



Formula/Name 



Refs 



amorphous olivine 
amorphous pyroxene 
crystalline forsterite 
crystalline enstatite 
amorphous carbon 



MgSi04 
Mgo.sFeo.sSiOa 
Mg2Si04 
MgSiOa 
AGAR 



Dorschner et al. (1995'! 
Dorschner et al. (1995) 
Servoin & Piriou (1973) 
Jaaer et al. (199|) 
Zubkoetal. (199&) 



St. Petersburg database'^ jHenninget al ."1999'. see Table|3}. 
Absorption efficiencies (Qats) are then calculated for spheri- 
cal, homogeneous grains of different sizes, representing com- 
pact grains, using Mie theory (.Mie,.1908i) . Interstellar dust, 
however, is not likely to be spherical and images of inter- 
planetary dust particles (see Bradley '2003", for review) sug- 
gest that nebular dust likely consists of porous and irregularly 
shaped aggregates. Additionally, the positions of the features 
in the mass absorption coefficients of crystalline silicates de- 
rived via Mie theory do not agree with observ ations of as- 
tronomical silicates (e.g., Bouwman et al. "2001!). Therefore, 
to simulate porous or irregularly shaped particles, we calcu- 
late opacities for continuous distributions of hollo w spheres 
(DHS), using the statistical method described in iMin et all 
(12005). In the Mie method, the particle can be represented by 
a homogeneous, filled sphere. In the DHS method, a distri- 
bution of hollow spherical silicate shells is used such that the 
vacuum volume filling fractions vary, but the masses remain 
the same as that of the homogeneous, filled sphere. The DHS 
method averages the scattering and absorption/emission cross 
sections of the set of hollow spheres and these cross sections 
can be compared to those of the homogeneous filled sphere 
of the same mass. The optical properties of individual large 
ho llow sph eres are calculated using the algorithm proposed 
by iToon & AckermanI (119811) for coated spheres, which was 
found to be extremely accurate (to 4 decimal places) when 
compared to measured efficiencies. We apply a simple analyt- 
ical application of the Rayleigh approximation as described in 
iMin et al. ( 2003 ) for small grains with x — 2'Ka/X << 1 and 
I ma: I << 1, where a is the radius of the particle and m is the 
refractive index. 

Figures I^Hll show the absorption efficiencies (Qabs) for 
grains of amorphous olivine, pyroxene, and crystalline 
forsterite for a range of sizes between 0.1 and 10 fim. Here 
we have used the grain radius (a) of the filled, homogeneous 
sphere as a proxy for the grain mass, which is held constant 
for both Mie and DHS models. As indicated in Figure|6l the 
slope of the overall spectrum varies significantly as a function 
of grain size. Thus, in order to compare the model Qabs to 
the normalized Spitzer spectra, we use a "continuum" fitting 
method which is similar to that applied to the Spitzer spectra. 
Although this method will most likely not find the physical 
continuum level, treating the observed and modeled spectra in 
the same manner will enable a more robust comparison. Nor- 
malized Qabs in the 10 and 20 fim regions for the modeled 
grain opacities of amorphous olivines calculated using Mie 
and DHS are shown in Figure |5j;,d and Figure |5^,f, respec- 
tively. The comparison of these normalized absorption effi- 
ciencies to the observed spectra involves an implicit assump- 
tion about the dominance of silicate emission in the 8-35 fim 

http://www.astro.uni-jena.de/Laboratory/Database/databases.html 
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Fig. 6. — Absorption efficiencies for amorphous olivine for grain sizes of 
0.1, 2, 5, and 10 fim. In the upper left panel, Qabs for o, = O.lfim has been 
multiplied by a factor of 5. The dashed lines show Qabs calculated using 
Mie theory for spherical grains. The solid lines show Qabs calculated using 
the distribution of hollow spheres (DHS) method described in the text. 
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Fig. 7. — Absorption efficiencies for amorphous pyroxene for grain sizes 
of 0.1, 2, 5, and 10 ^m. In the upper left panel, Q„hs for o, = O.l/^m has 
been multiplied by a factor of 5. Lines are as in Figurel6l 



region. Disk models (e.g.. lPollack et al ] ll994ft based on so- 
lar elemental abundances find that the main refractory com- 
ponents are olivine, pyroxene, carbonaceous material and/or 
organics, water ice, troilite (FeS) and metallic iron. They find 
that silicates clearly dominate the opacity for T > 500 K, but 
at lower temperatures the shapes and strengths of the 10 and 
20 fim features can be affected by the opacities of organics 
and water ice. The silicate emission studied here arises from 
the warm surface layers of these disks and ices are unlikely 
to play any role. Additionally, the excellent correspondence 
between the modeled silicate opacities and the observed spec- 
tra (cf. Figure|5|l, in both contrast and spectral shape, suggest 
that we are in a regime in which this comparison is valid. 

4.3. Grain growth 



C2D IRS I: Silicate emission and grain growth 



9 




10 20 30 10 20 30 



Fig. 8. — Absorption efficiencies for crystalline forsterite for grain sizes of 
0.1, 2, 5, and 10 fj,m. In the upper left panel, Q„,hx for i = 0-lA*m has been 
multiplied by a factor of 5. Lines are as in Figure|6| Note that the results of 
Mie/DHS theory for large (a > 2 fj.m) forsterite grains have not been con- 
firmed by laboratory studies. Additionally, emission feature positions are off- 
set for Mie and DHS models. DHS models have been found to more closely 
agre^ with astronomical observations of siUcates in disks tvan Boekel et aL 



Both the shape and strength of the silicate features are de- 
pendent upon grain size, with Mie and DHS models possess- 
ing slightly different silicate features. In this section we will 
compare the 10 and 20 fim silicate features in the observed 
spectra to the Mie and DHS models (Figure|5}- 

Both Mie and DHS models show that the 10 pim silicate 
feature strength decreases, and the feature flattens, as a func- 
tion of grain size. The 10 fim feature strength for Mie models 
(Figure |5j;) is largest for grain sizes of a « 0.1 /im. When 
a = 8 fim, the feature has almost disappeared. A similar trend 
can be seen in the models of hollow spheres using DHS meth- 
ods in Figure |5^. The shape of the 10 /im feature is slightly 
more round for the DHS than Mie models, however, and the 
peak strength is underestimated (with respect to Mie models) 
for smaller grain sizes (a = 0.1 /im) and overestimated for 
larger grain sizes (a — 6.0 /j,m). 

For the 20 /im feature (Figure |5ll,f), again we see that for 
both Mie and DHS models, the feature strength decreases and 
broadens when the grain size is increased from 2.5 to 5.0 /im. 
The emission feature also appears to shift to longer wave- 
length with increasing grain size, peaking near 18 /im (sim- 
ilar to our observations) for Mie and DHS models of amor- 
phous olivines with grain sizes of a = 3-5 /im. The shapes 
of the 20 fim features of the Mie and DHS models are also 
slightly different, with DHS models showing a much flatter 
feature. Additionally, the peak strengths of the DHS models 
are smaller than those of the Mie models for smaller grain 
sizes (a = 2.5 /im) and larger than the Mie models for larger 
grain sizes (a = 4.0 /im). 

The observed spectra (Figure |5t,b) are quite similar to the 
models, differing slightly for the weakest emission features. 
The weakest 10 /im features observed toward several sources 
are similar to the models of 5 /im amorphous olivine grains, 
but appear to be slightly narrower and more sharply peaked at 
the edges of the feature. Additional emission from crystalline 
silicates may be necessary to reproduce the "boxy" structure 



of such features. The strongest observed 20 /im features are 
similar in strength to models of amorphous olivine with grain 
sizes >2 /im. The weakest features are similar in strength 
to models of amorphous olivine with sizes less than approxi- 
mately 4.5 /im. 

Thus, the comparison of the observed spectra with mod- 
eled amorphous olivine opacities provides a qualitative un- 
derstanding of the effect of grain growth on the 10 and 20 /im 
features observed for our sample. In the following section, we 
will perform a more quantitative analysis. 

5. STATISTICAL ANALYSIS OF THE 10 AND 20 fjM FEATURES 

Previous observations of disks around HAEBE and T Tauri 
stars have shown that plots of the shape vs. strength of the 
10 /im emission feature may be indicative of grain growth. In 
§5.1 we interpret variations in the shape and strength of the 
observed silicate features in terms of source-to-source varia- 
tions in grain size. Because it is important to understand the 
uniqueness of this interpretation, in §5.2 we evaluate the de- 
pendence of the 10 and 20 /im feature strengths on several 
other grain properties, including the grain size distribution, 
relative sizes of grains emitting at 10 and 20 /im, grain com- 
position, crystallinity, and porosity. 

5.1. Strength-shape trends 

We find that the shapes of the continuum normalized 
10 /im emission feature from HAEBE and TTs are propor- 
tional to the feature s trengths ( as was previous noted by 
'van Boekel et al."2003', and later 'Przv godda et al.ll2()()3l and 
Kessler-Silacci et al. 2005). In Figure [9^, we plot the shape 
of the 10 /im feature versus the feature strength, S^'^j^J!^. As a 
proxy for the feature shape, we use the ratio of the normalized 
flux at the peak of the crystalline silicate feature relative to 
that of the amorphous silicate feature (S'n .a / S'g.g). The fluxes 
are integrated over regions of ±0. 1 /im around the central 
wavele ngth. Data fro m the li terature (lit, ^ an Boekel et 3l\ 
2003; Przvgodda et aL.2003: .Kessler-Silacci et al„2005 open 
symbols) are plotted in addition to the data collected in this 
study (c2d, filled symbols). Source parameters and original 
references for the lit sample are listed in Table |4] There are 
five sources that overlap between the c2d and lit samples: RU 
Lup, AS 205, DoAr 24E, HD 98922, and HD 163296. The 
10 /im feature strengths and shapes of these sources are con- 
sistent to within 11 ± 7%. In the statistical analysis, we use 
only our data for the duplicate sources, except in the case of 
AS 205 for which we will use the Przvgo dda et al.l (l2003l) 
spectra for each source in the binary. There is a strong cor- 
relation for the entire sample o f 10 /im features, which is con- 
sistent with that first noted bv Ivan Boekel et aTl j2003h . This 
trend is similar to changes in siUcate emission due to grain 
size variations illustrated in Figure |5] with flatter, weaker 
10 /im emission features, characteristic of larger grain sizes, 
appearing in the upper left of Figure|9^. 

In order to interpret this trend more quantitatively, we per- 
form the same analysis for the continuum normalized model 
opacities discussed in §4.2. We include grains of amorphous 
olivine, pyroxene, and mixtures of the two for grain sizes 
between 0.1 and 20 utn . The Maxwell-Garnett mixing rule 
jMaxwell GarnettI 1 1904ft is used to calculate effective cross 
sections for all mixtures. The model data are overlaid on the 
observational data in Figure [9^, with lines connecting points 
of varying grain size for each model. Note that the peak 
strength corresponds to different sizes depending on the dust 
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TABLE 4 

Stellar/Disk PARAMETERS for MIR spectra from the literature 





RA'^ 


DEC^ 


Spect ^ 


Age 


lIUliZi VV 


opcctrcil 




Source 


\J zuuu ) 


J zuuu ) 


Kei. 


(Myr) 


(A) 


Type 


Kei. 


GG Tau'' 


04 32 30.31 


+ 17 31 41.0 


Prz 


1.7 


31-54 


K6-K7 


1,2,3 


AA Tau 


04 34 55.45 


+24 28 53.7 


K-S 


0.9-2.4 


37.1 


K7 


1,2,4,5 


Lkca 15 


04 39 17.80 


+22 21 04.5 


K-S 


2-11 


23.1 


K5 


1,4,6 


DR Tau 


04 47 05.48 


+ 16 58 42.1 


Prz 


3.8 


95.4 


K4 


i,2 


GM Aur 


04 55 10.90 


+30 22 01.0 


K-S 


0.9-1.8 


96.5 


K7 


1,2,4 


SUAur 


04 55 59.38 


+30 34 01.5 


Prz 


3.0 


3.5 


G2 


2,7 


GWOri 


05 29 08.39 


+ 11 52 12.7 


Prz 




27-29 


G5 


2,3 


CR Cha 


10 59 06.97 


-77 01 40.3 


Prz 




30-44 


K2 


3,8 


TWHya 


11 01 51.91 


-34 42 17.0 


Prz 


9-10 


194 


K7 


1,3,9,10 


Glass l"" 


11 08 15.41 


-77 33 53.5 


Prz 


2.5 


22/1.4 


K4 


11-14 


WW Cha 


11 10 00.7 


-76 34 59 


Prz 


0.3-0.8 


66-67 


K5 


3,4,8,11 


Hen 3-600 A 


11 10 28.86 


-37 32 04.8 


K-S 


10 


12.5 


M4 


13,14 


IRAS 14050-4109 


14 08 10.3 


-41 23 53 


K-S 




2.0 


K5 


13,15 


RULup 


15 56 42.31 


-37 49 15.5 


Prz, here 


0.1-0.5 


136-216 


K7-M0 


3,8,16 


AS 205 NE 


16 11 31.40 


-18 38 24.5 


Prz, here 


0.1 


154.6 


K5 


2,17 


AS 205 SW 


16 11 31.40 


-18 38 24.5 


Prz,here 


<0.1 


54.9 


K5 


2,17 


DoAr 24E'' 


16 26 17.06 


-24 20 21.6 


Prz, here 


1.5 


5 


KO 


17,18 


Haro 1-16 


16 31 33.53 


-24 27 33.4 


Prz 




59-76 


K2-K3 


2,8,18 


AK Sco 


16 54 44.85 


-36 53 18.6 


Prz 


1.0 


3-9 


F5 


3,13,19 


SCrANW 


19 01 08.60 


-36 57 20.0 


Prz 


3.0 


73. 


K3 


13,17 


S CrA SE 


19 01 08.60 


-36 57 20.0 


Prz 


1.0 


61 


MO 


13,17 



MWC 480 


04 58 46.27 


+29 50 37.0 


K-S 


5-6 


18.3 


A2-A3 


1,20-23 


UXOri A 


05 04 29.99 


-03 47 14.3 


vB 


3-5 


20.0 


A2-A3 


1,19,20-23 


HD 37357 


05 37 47.08 


-06 42 30.2 


vB 






AO 


22 


HD 37806 


05 41 02.29 


-02 43 00.7 


vB 


6.3 




B9-A2 


22-24 


HD 95881 


11 01 57.62 


-71 30 48.4 


vB 




21.1 


A1-A2 


21,22 


HD 98922 


11 22 31.67 


-53 22 11.5 


vB,here 




27.9 


B9 


21,22 


HD 101412 


1 1 39 44.46 


-60 10 27.7 


vB 




14-20 


B9.5 


13,21,22 


HD 104237 


12 00 05.08 


-78 11 34.6 


vB 


2 


24.3 


A4 


13,20-23,25 


HD 142666 


15 56 40.02 


-22 01 40.0 


vB 


10 


0.8-3 


A7-A8 


13,20-22 


HD 144432 


16 06 57.96 


-27 43 09.8 


vB 




5-9 


A7-F0 


13,21-24 


HD 150193 


16 40 17.92 


-23 53 45.2 


vB 


3-5 


5 


A0-A4 


20-22,25 


HD 163296 


17 56 21.29 


-21 57 21.9 


vB,K-S,here 


5 


14.5 


A0-A2 


1,21-23,26 


HD 179218 


19 11 11.25 


+ 15 47 15.6 


K-S 


0.1 


18.2 


B9-A0 


21-23 


WW Vul 


19 25 58.75 


+21 12 31.3 


K-S 


>10 


18-24 


B9-A3 


3,4,20-22 


HD 184761 


19 34 58.97 


+27 13 31.2 


K-S 






A8 


27 



RE FERENCES. — (1) "Thi et al."2001> (2) 'Cohen & Kuh?'1979> (3) 'ReiDurth et al."199#i (4) 'Hartmann et al."1998t (5) 
Stometal. 1989; (6) Poncet et al. 1998; (7) Levreaiilt 1988; (8) Anpenzeller et al. 1983; (9) Webb et al. 1999; (10) Tones etall 
^m.; (11) Lawsonetal. 1996; (12) Gauvin & Strom 1992; (13) Greaorio-Hetem et al. 1992; (14) Chen et al. 1997; (15) 
Gregorio-Hetem & Hetem 2002; (16) Hughes et al. 1994; (17) Prato et al. 2003; (18) Bouvier & Appenzellei 1992; (19) 
Herbig & Bell 1988; (20) Nattaetal. 1997; (21) Acke et al. 2005; (22) The et al. 1994; (23) van den Ancker et al. 1998; (24) 
Meeusetall200lt (25)|van den Ancker et al. l99Tt (26)|Finkenzeller & Mundt 1984; (27) Miroshnichenko et al. 1 998 

Note. — Sources in the top portion of the table are T Tauri stars; sources in the bottom portion are HAEBE stars. 

'^The RA and DEC quoted here were obtained from SIMBAD and do not necessarily represent the position of the referenced MIR 
spectroscopic observations. 

^Re ferences for silicate spectroscopy: vB — Ivan Boekel et aH2003i Prz — IPrzygodda et aH2003t K-S — IKessler-Silacci et alj 
I2005t here — this paper. 

"^All Ha lines are in emission 

'^This source is a binary that is unresolved in the referenced MIR spectroscopic observations. Stellar/disk parameters quoted from 
the literature include both sources. 



composition and whether DHS or Mie models are used (see 
bottom panel of Figure|9^). 

The observed trend is best fit with models of amorphous py- 
roxene/olivine mixtures. Models of pure amorphous olivines 
calculated using DHS are inconsistent with the observations. 
These models possess generally larger Sn.a-to-S'g.s ratios 
than the data, particularly for strong features. Mie mod- 
els of pure amorphous olivines fit the data much better, but 
begin to deviate from observations near 5^^^™ « 2. The 
slope of the observed trend can be best matched using mix- 
tures with 30% amorphous olivines and 70% amorphous py- 
roxenes. This pyroxene fraction is much larger than that in- 



ferred for diffuse ISM grains in the galactic center (15.1% 
pyroxene; Kemper et al. 2004), indicating a substantial con- 
version of olivines to pyroxenes in young stellar environ- 
ments. Enhanced pyroxene-to-olivine abundances were pre- 
viously observed toward high-mass protostars with ISO and 
may be explained by He+ sputtering in high-velocity shocks 
(JDemvk et al. 2001, and references therein). 

Models of olivine-pyroxene mixtures fit the data best for 
homogeneous filled spheres (Mie) and hollow spheres (DHS), 
but do not fit equally well for all grain sizes (peak strengths). 
For the weakest 10 ^im features (5"^"^™ < 1.75), hollow 
sphere models indicate grain sizes of up to 1.3 times those of 
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Shape and strength of the a) 10 and b) 20 fim silicate features. The shapes of the features (probed by the normalized flux ratios, Sn.a/Sg 



and <S'23.8/'S'i9) are plotted versus the feature strengths (S^^'^^ 



and 5" 



-'peak ^ ''^^ '"P P™^''' °f each figure. The data observed in this study are denoted 
by filled squares (TTs) and filled triangles (HAEBE). Data for 10 fim silicate features from the literature ( van Boekel et al. 2003; Przvgodda et al. 2003| 
iKessler-Silacci et al. 2005) are always represented by open symbols. The black line in a) indicates the correlation found for the 10 fim features including 
all data, and follows the relation, y = A + Bx, where A = 1.37 ± 0.05 and B = —0.24 ± 0.03, with a correlation coefficient of r = —0.7. The 20 /im 
features also show a correlation (r = — O.5.), with A = 1.23 it 0.08 and B = —0.21 ± 0.06. Models of amorphous olivines, pyroxenes and mixtures of the 
two for various grain shapes are overlaid as dashed-dotted lines connecting points representing grains with a linear size distribution between 0. 1 and 20 ^im (see 
text for details). The bottom panels of each figure show the grain size vs. peak strength for the models. 



homogeneous spheres for features of the same strength. Both 
mixture models indicate that the largest changes in the shape 
of the feature occur over a small range of grain sizes, from 
approximately 1-3 fjm, with exact values depending on the 
model used. For the weakest silicate features (S'pg^™ < 1.7), 
large changes in grain size result in only small variations in 
the strength/shape of the silicate feature. Therefore, it be- 
comes increasingly difficult to determine exact grain sizes for 
the weakest silicate features. 

The data discussed above cannot be explained purely by 
size variations, and in some cases can indeed be attributed 
to crystalline silicates and/or PAH emission. For the largest 
grain sizes. Figure shows several spectra with 11.3-to- 
9.8 /im ratios that consistently lie above the modeled ratios. 
These data may indicate the presence of crystalline forsterite 
or PAH emission in the 1 1.2-1 1.3 /im region. For a cutoff at 
8113/8^,% = 1.0, we find that 13 sources in our sample sat- 
isfy this criterion, and thus require emission near 11.3 /xm that 
is in excess above models of amorphous olivines/pyroxenes, 
indicating that the features arise from crystalline silicates or 
PAHs. These sources are marked with a "Y" in the column 
3 of Table El Of these 13 sources, 4 show clear evidence of 
crystalline silicate emission at A > 25 /im (DoAr24E, GY 23, 
HD 98922 and DL Cha). Clear evidence for PAH emission 
features in the SL spectra can be seen in 3 additional sources 
(LkHa 330, IRAS 03446H-3254, Haro 1-17 and HD 135344). 
The spectra of 3 more sources (VSSGl, LkHa 327, and RR 
Tau) show evidence of both PAH and crystalline silicate fea- 
tures. The remaining 3 sources with Sn.a/Sg.g > 1-0 (T 
Cha, SR 21, and HD 135344) display narrow isolated 1 1.3 /im 
emission and have very simple spectra that show only shal- 
low ^^20 /im silicate emission. This likely indicates that the 
11.3 /im features in these 3 spectra arise from PAH, and not 
crystalline silicate, emission. Therefore, for all of the spec- 
tra with 5*11. 3/5*9. 8 > 1-0, the additional flux at 11.3 /im can 
indeed be explained by crystalline silicate and/or PAH emis- 



sion. The nature of the 11.3 /im features and relationship to 
PAH emission will be discussed more thoroughly in Geers et 
al. (in prep.). Aside from the spectra with S'n.a/S'g.s > 1-0, 
the data presented in Figure |9^ are consistent with grain size 
variations. 

The 20 /im features show a shape-strength trend which is 
very similar to that seen for the 10 /im features (Figure|9j)). 
Here we use the ratio of the flux at 23.8 to 19 /im, again the 
peaks of the most prominent amorphous and crystalline sili- 
cate features, as a proxy for the feature shape. This plot in- 
cludes only c2d data, as continuum normalized ISO data were 
not available. The paucity of spectra with strong 20 /im fea- 
tures may be related to the prominence of weak, flat 10 /im 
features in the c2d sample, indicating large grain sizes. Al- 
though the sample size (and range in peak strength) is smaller 
for the 20 /im features, the shape-strength trend is still promi- 
nent. This trend possesses a similar slope and y-intercept to 
that seen for the 10 /im features and is again most closely 
matched with the Mie model of an olivine/pyroxene mixture, 
which agrees to within S'23.8/'S'i9 = 0.01, or ^1%. There 
appears to be little difference in grain composition between 
grains being probed by the 10 and 20 /im features. 

All of the 20 /im features are consistent with models of 
amorphous olivine/pyroxene mixtures of sizes between 0.1 
and 20 /im, with no deviations indicative of strong crystalline 
features. A similar fit is obtained for the Mie and DHS mod- 
els of the olivine/pyroxene mixture, but for the DHS mod- 
els much larger grain sizes (a > 20 /im) are needed. The 
data are also consistent with DHS and Mie models of amor- 
phous olivines (to within S'23.8/5'19 = 0.08, or ~8%), for 
much smaller grain sizes (a<5 /im). This makes exact grain- 
size determination difficult. Furthermore, a large scatter in the 
modeled data points can be seen for S'p"^™ < 1 . 1 due to dif- 
ficulties in fitting the continuum for such weak features. This 
prevents identification of small deviations in feature shape of 
between the models and data for large grain sizes (as were 
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TABLE 5 
Model parameters for FiouREfTol 



Panel 


Parameter 


Vai'iables 


Range 


b 


size distribution" 


P 


0.01-3.0 Aim 
10-100 fim 
4.5-2.5 


c 


sizes for 10 vs 20 /im 


P 

a-'"- 


4.5-2.5 
0.1,1.0 fim 
(0. 1-20) X 0.1 fim 
(0.5-3) X 1.0 Atm 


d 


composition 


a^" "="^^" 

rmn rmji 

olivine: carbon 
olivine:pyroxene 


0.01-3.0 fj.m 
100:0-50:50 
100:0-0:100 


e 


crystallinity 


a^" - a^" 

mm mm 

olivine:forsterite 
olivine:enstatite 


0.01-3.0 fj.m 
100:0-0:100 
100:0-0:100 


f 


porosity 


min mm 

olivine:vacuum 


0.01-3.0 iim 
100:0-20:80 



"Using a differential grain size distribution of dn{a)oca ^ da where the absorp- 
tion efticiency is defined as < Qahs >— j^^^^ Qabsdn{a)/ J^^'^^ dn{a). 



seen for the 10 /im features). However, larger deviations from 
the models are not seen. Thus, although there is evidence in 
some individual spectra for crystalline silicate features, the 
entire trend is consistent with amorphous olivine/pyroxene 
mixtures with source-to-source size variations. 

For the 48 sources observed in this study, the 10 and 20 fim 
features both exhibit strength-shape trends consistent with 
source-to-source grain-size variations, with most of the ob- 
served features similar to models of grains with sizes much 
larger than that of the ISM (a >> 0.1/im). This indicates 
that grain growth must occur quickly in these disks. Further- 
more, as the observed features arise from silicate emission 
in the disk surface layers, and the largest grains will gravita- 
tionally settle to the disk midplane, the emitting grains likely 
represent the low-size tail of the silicate grain size distribution 
within the disk. This effect is enhanced by the fact that larger 
grains have weaker silicate emission features. 

5.2. 10 vs. 20 feature strengths 
One of the most interesting aspects of this sample is a set 
of sources with very strong Sl^„[!'Y^ but weak S'^^'^^ . These 
"outliers" can be clearly seen in a plot of the strengths of the 
10 and 20 /im features in the observed spectra (Figure [Tol 
Table ??). The feature strengths appear to be correlated over 
the entire sample (r = 0.36, 96% probability of correlation). 
However, examination of Figure [lOk reveals that the bulk of 
the data is clustered around 5*^°^™ « 1.3 and 5*^°^™ « 1.5 
and ~11 outliers have 10 fim and 20 iim features with very 
different strengths. Furthermore, most of these outliers have 
strong 10 fim features (5*^°^^ > 2) and weak 20 /im features 
(^peak^ < 1.6). No sources in our sample have strong 10 fim 
features and equally strong 20 fim features. One may expect 
weak 20 fim features to be more abundant as most sources 
show weak 10 fim features (cf. Figure|9t) and grain growth is 
occuring quickly in these disks. However, the fact that there 
are no sources with strong 20 fim features and weak 10 fim 
features is significant. 

In order to interpret the distribution of peak strengths in 
the observed spectra, we examine the effect of a variety of 
grain parameters on the strengths of the 10 and 20 fim fea- 



tures. Factors affecting silicate feature strengths include, but 
are not limited to, the grain size distribution, relative sizes 
of grains emitting at 10 and 20 fim, grain composition, crys- 
tallinity, and porosity. We attempt to assess the i nflue nce 
of each of these parameters in panels b-f of Figure ^] and 
§5.2.1-5.2.4 for models of normahzed Qats calculated as de- 
scribed in §4.2. 

5.2.1. Grain size distributions 

The effects of grain size are explored by modeling opaci- 
ties for a set of grains with a differential grain size distribu- 
tion dn{a) oc a^P da, and varying the minimum grain size 
amin, maximum grain size Umax, and power-law index p, as 
shown in Figure 1 10b . All three parameters have very simi- 
lar effects; increasing amin or amax or decreasing p results 
in more large grains and reduces both the 10 and 20 fim fea- 
ture strengths, but does not affect the ratio of the two. Thus, 
variation of the grain size distribution cannot account for the 
extreme lO-to-20 fim ratios of the outliers in Figure [Tok. The 
bulk of the observed data set is centered around 10 and 20 fim 
feature strengths similar to that of a grain size distribution 
with amin = 1.0 fim, amax = 100 fim, and power-law in- 
dex p = 3.5. 

The emission at 10 and 20 fim likely comes from differ- 
ent populations of grains. This is a simple consequence of 
the temperature dependence of the emission as a function of 
wavelength (e.g., emission at 20 fim can arise from colder 
grains than emission at 10 fim). We therefore expect the con- 
tribution of a particular grain to the 10 and 20 fim features 
to depend on the location of the grain. In order to deter- 
mine the relative contribution of grains from different radii 
to the 10 and 20 fim features, we perform an exercise using 
a sim ple two-layer rnodel (CGPLUS. Dullemond et al. 200 U 
iChiang & GoldreichllT997i) to calculate the cumulative disk 
flux as a function of radius at 10 and 20 /im for a typical 
T Tauri star. Figure shows that the radii probed by the 10 
and 20 fim silicate features can indeed be quite different, with 
the 20 fim feature probing radii up to 10 times those probed 
by the 10 fim feature. 

In Figure llOt . the relationship between the sizes of the 
grains emitting at 10 and 20 fim is explored. We vary amin 
for the models of the 20 fim feature and keeping amax, P, and 
cimin for the models of the 10 fim feature constant. The entire 
distribution of the data can be reproduced by these models. 
Most data correspond to a'^'^™ = 1 fim and a^°^ = 0.1- 
3 fim, but the data in the right side of Figure llOt require 
large differences between the sizes of grains emitting at 10 
and 20 fim (a^°f " = 0.1 fim, al^^l" = 2-3 fim). 

The 10 and 20 fim features are clearly probing different dust 
populations. The bulk of the data can be reproduced by mod- 
els of 1.0 fim grains, indicating significant growth in these 
disks compared to the ISM. Although the data are clustered 
centered around the dashed lines that denote equal grain sizes 
for 10 and 20 fim emission, many spectra can be better repre- 
sented by much larger grains emitting at 20 fim than at 10 fim. 
This may be explained, if the emission is optically thin, by a 
difference in the t — 1 surface depth as a function of wave- 
length of the emission; 20 fim features probe a deeper layer of 
the disk where grain sizes are likely larger due to larger den- 
sity and dust settling. One would also expect, however, that 
dust at larger radii (and lower temperature), which is likely 
smaller due to decreased densities, contributes more to the 
20 fim feature than the 10 fim feature. There is thus a com- 
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Fig. 10. — Relative strengths of the 20 versus 10 fim features. Models are oveiplotted showing variations in b) grain size distribution, c) ratios of grain sizes 
emitting at 10 and 20 ^m, d) grain composition, e) crystallinity, and /) porosity. The reference opacities are calculated via Mie theory for pure olivine grains 
with differential grain size distributions with a„ii„ = l.Ofim, amax = 100/^m, and power-law index p = —3.5 for both the 10 and 20 ^m features. Arrows 
in panel b indicate the direction of increasing slope p, a^in and imax- Gray shaded regions in c) indicate the range in 10 and 20 fim feature strengths covered 



by variations in a^"^™ 



with a 



■ 0.1 um and a^'''"" 



1.0 ^m. Arrows in panel c indicate the direction of increasing lO-to-20 ^m grain size ratios over 



the gray shaded regions. Gray shaded regions in d-f indicate the range in 10 and 20 fim feature strengths covered by variations in x/olivine ratios, where x is d) 
pyroxene or amoiphous carbon, e) forsterite or enstatite, and /) vacuum for grains of constant size. Arrows in d-f indicate the direction of increasing x/olivine 
ratios. Symbols are defined as in Figure l9l 
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Fig. 1 1. — Cumulative flux at 10 ^m (solid line) and at 20 iim (dashed 
line) as a function of radius. Fluxes are calculated using the CGPLUS model 
jpullemond et al. 2001) for a disk around a T Tauri star with T* = 4000 K, 
L* = 0.58 Lq containing olivine grains of a differential grain size distri- 
bution with amin = 1 Mrn, cimax = 100 /im and p = —3.5. 95% of the 
10 /im flux comes from within ~1 AU, while 95% of the 20 fim flux comes 
from within ~10 AU. 



peting effect between the disk height and radius being probed. 
This may explain why most of the data points are centered 
around the dashed lines that denote equal grain sizes for 10 
and 20 /im emission. In the case of the outliers, enhanced 
dust settling may make the disk height differential between 
dust probed by the two features more important than the dif- 
ference in radius. 

5.2.2. Grain composition and crystallinity 

Grain composition and crystallinity also affect the strengths 
of the silicate emission features, resulting in ranges of peak 
strengths that are very similar to those produced by varying 
the grain sizes (with 0,"^'^™ — a^j>'")- Thus, effects of grain- 
size variation, composition and crystallinity are difficult to 
disentangle using lO-to-20 /im flux ratios. 

In Figure [Toh. we vary the percentage of amorphous py- 
roxene and amorphous carbon relative to amorphous olivine. 
To simplify matters, we use the same composition and grain 
sizes for the grains emitting at 10 and 20 ^m. Increasing the 
pyroxene-to-olivine fraction has the largest effects for small 
grains (<1.0 /im), resulting in an increased 10 /im feature 
strength and roughly constant 20 fim feature strength (as seen 
previously in Figure|9jl. Increasing the carbon-to-olivine frac- 
tion has an even larger effect on grains of small sizes, with a 
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larger effect on the 20 /im feature than the 10 /im feature, as 
the peak-to-continuum ratio is decreased due to an increase in 
the effective "continuum." 

Variations in crystallinity (using DHS models) also result in 
the same range of peak strengths, but highly crystalline grains 
can produce features with the same peak strengths at smaller 
grain sizes ("Figure [Tofe). When the fraction of forsterite or 
enstatite is less than ^50%, the feature strengths remain very 
similar to amorphous olivine. Beyond a crystalline fraction 
of ^^50%, the strong and narrow crystalline emission features 
begin to dominate and the strengths of the peak emission in 
the 10 and 20 fim regions appear to increase dramatically 
for the same grain size distribution. For models of >50% 
crystalline forsterite and enstatite, the gray shading in Fig- 
ure llOb extends far beyond the plotted range of the y-axis 
(^peak^ >> 2.28). Thus, for primarily crystalline grains, 
grain size distributions with amin < 2.0 fxm can account for 
the entire gray area shown. 

Differing composition and crystallinity can account for the 
10 and 20 /im ratios of the bulk of the observed spectra, but 
not the outliers. Additionally, the 10 and 20 /im ratios of mod- 
els of varying grain composition and crystallinity encompass 
ranges in which no observational data lie. In particular, spec- 
tra showing weak 10 fim features (S^'^^J!^ > 1.7) and strong 
20 /im features can be easily produced with models of small- 
sized {amin < 1-0 /im) amorphous silicate/carbonaceous 
grains or moderately-sized {a„iin < 3.0 /im) crystalline sili- 
cate grains, but they are not seen in the c2d sample. These ab- 
sences confirm the idea that the 20 /im feature probes regions 
of larger grain size than does the 10 /xm feature, irrespective 
of grain composition and crystallinity. 

5.2.3. Dust porosity 

Finally, in Figure we explore the effects of grain opac- 
ity by increasing the vacuum fraction of a set of olivine 
grains and by using the distribution of hollow spheres method 
(DHS). In both cases, we integrate over a grain size distribu- 
tion as in Figure [Toh-e. The primary difference between the 
two methods is that the opacities calculated with DHS are the 
average of a set of hollow spheres of varying vacuum frac- 
tions for a given mass, which is translated to grain size. The 
Qabs are then integrated over a grain size distribution. For 
the olivine-Hvacuum grains, the opacities and Qabs are calcu- 
lated and integrated over grain size for each volume fraction 
individually. The variation over vacuum fraction (light gray 
in FigureFTOI^ does not have a large effect on the lO-to-20 /im 
ratio and is consistent with the DHS 10 and 20 /im strengths 
and variation with amin- Variations in dust porosity can only 
explain a fraction of the bulk of the sample, which overlaps 
with the regions covered by models of varying grain composi- 
tion, and cannot explain the outliers. Thus it is not necessary 
to consider porosity to explain the observed 10 and 20 /im 
feature strengths. 

6. GRAIN GROWTH AND STELLAR/DISK PROPERTIES 

In order to evaluate the dependence of the strength-shape 
trend noted in §5 on stellar and disk parameters, we plot the 
10 fim shape versus strength again in Figures [T2k-c. color- 
coded by spectral type, stellar age, and Ha equivalent width 
(for TTs only). Stellar/disk parameters for each source ob- 
served in this study {c2d) were collected from the literature 
and are shown in Table [D Stellar/disk parameters for 10 
data obtained from Ivan Boekel et al.l(l200l ■ iPrzygodda et al.l 



(120031), and lKessler-Silacci et alJ ( 120051) (lit) are shown in Ta- 
ble@] We evaluate the probabiUty that the sets of spectral type, 
stellar age or Ha equivalent width bins are drawn from the 
same distribution on the 10 /im feature shape-strength relation 
shown in Figure|9t by evaluating the two-sided Kolmogorov- 
Smirnov statistic (D)''* for each stellar/disk parameter As 
we are evaluating the dependence of the strength-shape trend 
on stellar/disk parameters, we use only the sources with 
"5*11. 3/5*9. 8 < 1-0, thus removing sources with substantial 
emission from crystalline silicates and/or PAHs. Some stel- 
lar/disk parameters are not available in the literature, and 
therefore not all sources are included in Figures I12h . b, and 
c. 



6.1. Spectral Type 

As the data in our sample (filled symbols) appear to be con- 
sistent with the largest silicate grain sizes, and probe disks 
around primarily K and M stars, we first looked for connec- 
tions between grain size and spectral type. We divide the sam- 
ple by spectral type into four groups: 1) A/B, 2) F/G, 3) K and 
4) M. The A and B stars and F and G stars were combined 
into single groups as few stars of each type were observed. 
The Kolmogorov-Smirnov (K-S) test is used to compare two 
arrays of data values and the process is repeated for each com- 
bination of groups 1 through 4. In general, there appears to 
be a dependence on spectral type, with A/B stars being clus- 
tered toward the middle and lower right of the 10 /im feature 
strength-shape trend and M stars being more clustered toward 
the upper left of FigurJT2l This can be more easily seen in 
Figure where histograms of the 10 /im shape-to-strength 
ratio for A/B stars and M stars are compared. The population 
of A/B stars peaks at ratios near 0. 1 and decreases for larger 
ratios, while only a few M stars have shape-to-strength ratios 
near 0. 1 and the population increases for larger ratios to peak 
near 0.5-0.6. 

Although it appears that 10 /tm silicate emission may be 
different in A/B stars versus M stars, a larger data set is 
needed to confirm this result. The calculated probabilities 
that these groups are drawn from the same distribution range 
from 15% (D = 0.29) for spectral types A/B and M to 98% 
(D = 0.17) for specti-al types A/B and F/G. Objects that ai-e 
close in spectral type, in subsequent groups (e.g., 1 vs. 2, 
2. vs. 3), have large probabilities (69-98%) of being drawn 
from the same distribution. Groups with larger differences in 
spectral type (e.g., 1 vs. 3, 2. vs. 4) have smaller proba- 
bilities (15-50%). Of all these comparisons, however, only 
the K-S test for group 1 vs. 4 (A/B vs. M) indicates a low 
enough probability (< 20%. lPress et alJl992 ) to suggest that 
the differences in the populations are statistically significant. 
Furthermore, the number of objects in each group is still quite 
low (13 A/B stars, 4 F/G stars, 24 K stars, and 15 M stars) 
and the spectral types typically have errors on the order of 3 
subclasses. 

The dependence of the strength-shape trend with spectral 
type may be explained in terms of the disk temperature and 
density structure. If we assume that the 10 /im feature always 
probes grains of roughly the same temperature, then the radius 
being probed by the 10 /im emission will change as a function 
of the temperature structure of the disk. The disk temperature 

The two-sided Kolmogorov-Smimov statistic is the difference between 
the cumulative distribution function (CDF) of t wo sets of data a nd is calcu- 
lated using KSTWO (from Numerical ReciDes: |Press et aH1992ft . 
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Fig. 12. — Shape and strength of the 10 silicate features. The shape of the feature (S11.3/ Sg.s) is plotted versus the feature strength (5^1^^™) in all panels. 
The points are color-coded by a) spectral type, b) age and c) Ha equivalent width. Silicate features from the literature (van Boekel et al. 2003; Przvaodda et al] 
l2003i: .Kessler-S ilacci et al. 2005 ) are always represented by open symbols. The black line in all panels indicates the congelation found for the 10 /^m features 
shown in FigureWt including all observations. This trend is not strongly dependent on HqEW or age, but may be related to spectral type. 
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Fig. 13. — Histogram of the 10 ^m feature shape-to-strength ratios for the 
A/B and M stars shown in Figure fT2t . Spectra of A and B stars show mor e 
centrally peaked 10 fim features, lying in the middle and left of Figure [T2t . 
while spectra of M stars show primarily flattened features, lying in the upper 

left of Figure ll2t . The bin widths are 0.1 in units of (5ii.3 /Sq.s) I (S^J^^). 



T{R) varies as a function of radius as 



T{R) 



(from'Dull emond & Dom inik 2005). where a is the disk flar- 
ing angle and Ri, and are the stellar radius and temperature. 
R is the radius corresponding to temperature T{R), in this 
case the dominant radius being probed by the silicate emission 
feature, which we will call Rsu- Assuming that T{Rsii) and a 
are constant over the entire disk sample, and using typical val- 
ues of i?* and for A/B and M stars, Rsu^a/b I Rsii,M~ 10- 
25. So the radius probed by the 10 /im emission feature is 
quite different for stars of spectral types A/B and M. Further- 
more, the density generally increases with decreasing radius 
and grain growth increases at higher density, while dust set- 
tling decreases at higher density. Thus, the combination of 
faster grain growth and slower dust settling at smaller radii 



could lead to the larger grains that we see in M stars versus 
A/B stars. Although this simple theory is consistent with the 
observed correlation of grain size with spectral type, several 
other spectral-type dependent factors (e.g.. X-rays, UV radi- 
ation, stellar/disk winds, etc.) likely influence grain size in 
these disks. 

6.2. Age 

Although one might expect grain growth and crystalliza- 
tion to occur over time, as the disk matures, we see no clear 
relation ship between the strength-shape trend and stellar age 
(Figure[2p). The ages for our sample are culled from the lit- 
erature and range from 0.5 to 6 Myr with errors of up to a few 
Myr (i.e., Hillenbrand & White 2004). We therefore divide 
up the sample into three approximately equally sized groups: 
stars with ages of 1) < 1.5 Myr, 2) 1.5^ Myr, and 3) >4Myr, 
and follow a similar method as described above for the analy- 
sis of spectral type. The strength-shape trend does not appear 
to be strongly related to the age of the star (Figure[T2b). with 
probabilities of >47% that all three groups are drawn from 
the same distribution. However, this analysis would be more 
conclusive for stars within one cluster with a well defined age, 
for which relative ages could be more accurately determined. 
With the available ages, these results indicate that the stel- 
lar age of PMS stars is not directly related to the state of the 
small, /im-sized grains in these disks. 

In addition to the difficulties in determining the disk age, 
the lack of an age dependence of the shape-strength trend 
can be understood when grain fragmentation and disk tur- 
bulence are considered. Dullemond & Dominik (.200 5 ) find 
that in models including only basic coagulation mechanisms 
grain growth occurs very quickly, such that the SEDs show 
much weaker IR excess (even at young ages of 1 Myr) than 
do disks around most classical TTs. The settling of small 
grains, removing them from the disk surface, occurs over 
even faster time scales. The lack of sources with strongly 
peaked 10 /im features (and therefore small grains) in the 
substantial sample observed here is consistent with this fast 
timescale for grain growth (to 1 /im sizes). The presence of 
strongly peaked 10 /im features in some sources (particularly 
in the lit sample) may be indicative of the replenishment of 
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small grains in disk surface layers through mechanisms such 
as fragmentation and/or turbulent mixing. (Some amount of 
large grains will also be mixed into disk surface layers). The 
models presented in Dullemond & Dominik (2005) suggest 
that such mechanisms may allow disk systems to maintain 
quasi-stationary grain size distributions, with grain growth 
and small-grain replenishment rates in equilibrium for long 
periods of time (>1 Myr). Thus the typical grain sizes found 
in disks around pre-main sequence stars are likely not to be di- 
rectly related to disk age, but indicative of a number of other 
factors, such as the strength of the turbulence and the gas mass 
in the disk. 

6.3. Accretion 

We also do not see a correlation of the strength-shape trend 
with the equivalent width of the Ha emission lines (HaEW), 
probing disk evolution (Figure 1121). The Ha equivalent 
width is a tracer of active stellar accretion from the inner disk. 
Thus clearing of the inner disk will result in lower HaEW. 
In Figure I12t . we separate the TTs in our sample into two 
groups using the equivalent width of Ha: 35 classical T Tauri 
stars (cTTs; HaEW > 10 A) and 14 weak-lined T Tauri stars 
(wTTs; HaEW < 10 A). Although there is no clear observa- 
tional evidence that cTTs are progenitors to wTTs, wTTs are 
often described as cTTs in which the inner disk has dissipated. 
Therefore, we might expect that the clearing of the inner disk, 
and thus the transition between the classical and weak-line 
T Tauri star phases, would be related to grain growth and 
the 10 fim feature strength-shape trend. We find that there 
is a high probability (^70%) that the two groups are drawn 
from the same distribution, however, indicating that there is 
no clear relationship between HaEW and grain growth in 
these disks. This indicates that either HaEW is not a good 
probe of the disk evolutionary state or stellar accretion rate or 
that grain growth is not related to these quantities. We also 
note that Ha emission can be quite variable and the values 
used here are time averages, using the full range of observed 
HaEW as error bars. The relation between HaEW and sili- 
cate feature shape-strength would more precisely tested with 
simultaneous MIR and optical spectroscopy. 

7. CONCLUSIONS 

Spectra in the ^^5-35 fim region have been obtained for 
disks around 40 TTs and 7 HAEBE stars using the Spitzer 
Space Telescope, as a subset of the c2d IRS survey. This is 
the first significant sample of T Tauri stars for which both 10 
and 20 /im features are available. Broad silicate features at 10 
and 20 fim, similar to emission from amorphous silicates, are 
prominent toward the observed sample of TTs and HAEBE 
stars. Emission features from crystalline silicates are also evi- 
dent in the observed spectra, with the most prominent features 
visible in the 33-36 region. 

We performed a statistical analysis of the shapes and 
strengths of the 10 and 20 /Ltm features and find the follow- 
ing: 

— If the full IRS spectra are considered, the data are most 
consistent with source-to-source variations in grain 
size, with the bulk of the sources in our sample indi- 
cating sizes of 1 fim or greater. This and the lack of 
strongly peaked 10 /im sources in our sample of 40 



sources is consistent with fast grain growth (from 0.1- 
1.0 fim sizes) in the surfaces of these disks. 

— Long wavelength (33-36 /im) crystalline silicate fea- 
tures are seen toward ^half of the TTs disks in our 
sample. This indicates that significant dust processing 
is also occuring in these disks. Only 13 of these spec- 
tra also exhibit 1 1.3-to-9.8 /im ratios that cannot be re- 
produced by models of amorphous silicates, suggesting 
that the degree of crystallinity deduced from the 10 //m 
region alone is underestimated. 

— A subset of the observed spectra have particularly large 
lO-to-20 /im ratios, which can only be reproduced 
with much smaller grains emitting at 10 /im (a„ji„ = 
0.1 /im) than at 20 /im {amin = 2 — 3 /im). This can 
be explained if the 20 /im emission arises from deeper 
in the disk than the 10 /im emission in disks where dust 
settUng has occurred. 

— The 10 /im feature strength vs. shape trend is not corre- 
lated with age or disk evolutionary state (HaEW). This 
suggests the importance of turbulence and the regener- 
ation of small (/im-sized) grains on the disk surface. 

— The 10 /im feature strength vs. shape trend is related to 
spectral type, with M stars showing significantly flatter 
silicate features (larger grain sizes) than those of A/B 
stars. This may be related to a difference in the ra- 
dius probed by the emission, which should increase as a 
function of the disk temperature, and thus stellar lumi- 
nosity. However, the observed correlation could also be 
indicative of other spectral-type dependent factors (e.g. 
X-rays, UV radiation, stellar/disk winds, etc.). 

As the observed features arise from silicate emission in the 
disk surface layers, and the largest grains will gravitationally 
settle to the disk midplane, the emitting grains likely represent 
the small-size tail of the silicate grain size distribution within 
the disk. This simple picture of grain growth and settling is 
complicated by vertical mixing, which can bring both large 
and small grains back to the disk surface, and fragmentation, 
which results in the replenishment of small grains through- 
out the disk. Comparisons with probes of grain size in the 
disk midplane and inner disk clearing, as well as probes of 
disk turbulence and gas-dust interactions, are vital to under- 
standing the results presented here in the context of planet 
formation. 
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